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Abstract. It is shown that the lens equation for a binary gravitational lens being a set of two coupled real fifth- 
order algebraic equations (equivalent to a single complex equation of the same order) can be reduced to a single 
real fifth-order algebraic equation, which provides a much simpler way to study lensing by binary objects. 
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1. Introduction 

The gravitational lensing due to a binary system has at- 
tracted a lot of interest since the pioneering work by 
Schneider and Weifi (1986). The lens equation adopted 
until now is a set of two coupled real fifth-order equa- 
tions, equivalently a complex fifth-order equation (Witt 
1990) which is based on a complex notation introduced 
by Bourassa, Kantowski and Norton (1973, 1975). 

The number of images is classified by curves called 
caustics, on which the Jacobian of the lens mapping van- 
ishes on a source plane. Caustics for two-point masses are 
investigated in detail and locations of caustics are clari- 
fied based on a set of two coupled real fifth-order equa- 
tions as an application of catastrophe theory (Erdl and 
Schneider 1993) and on a complex formalism (Witt and 
Petters 1993), which is developed as an efficient method to 
compute microlensed light curves for point sources (Witt 
1993): In the binary lensing, three images appear for a 
source outside the caustic, while five images are caused for 
a source inside the caustic. For a symmetric binary with 
two equal masses, the lens equation for a source on the 
symmetry axes of the binary becomes so simple that we 
can find the analytic solutions (Schneider and Weifi 1986). 
In star-planet systems, the mass ratio of the binary is so 
small that we can find approximate solutions in general 
(Bozza 1999, Asada 2002). The approximate solutions are 
used to study the shift of the photocenter position by the 
astrometric microlensing (Asada 2002). 

Nevertheless, it is quite difficult to solve these equa- 
tions, since there are no well-established methods for solv- 
ing coupled nonlinear equations numerically with suffi- 
cient accuracy (Press et al. 1988). We show that the lens 



equation can be reduced to a single master equation which 
is fifth-order in a real variable with real coefficients. As a 
consequence, the new formalism provides a much simpler 
way to study the binary lensing. 

2. Lens Equation for a Binary System 

We consider a binary system of two bodies with mass Mi 
and M2 and separation vector L from the object 1 to 2. 
For a later convenience, let us define the Einstein ring 
radius angle as 
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where G is the gravitational constant, M is the total mass 
Mi + M 2 , and Dj^, Dg and D-^g denote distances between 
the observer and the lens, between the observer and the 
source, and between the lens and the source, respectively. 
We choose the position of the object 1 as the coordinate 
center. In the unit of the Einstein ring radius angle, the 
lens equation reads 



= 0- i/i 
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where (3 and 6 denote the vectors for the position of the 
source and image, respectively and we defined the mass 
ratio and the angular separation vector as 
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where we defined 



We have an identity v\ + = 1. For brevity's sake, v-i is 
denoted by v. Equation (Q) is a set of two coupled real 
fifth-order equations for (9 x ,9 y ), equivalent to a single 
complex fifth-order equation for 8 X + iO y (e.g. Witt 1990, 
1993). 

Let us introduce polar coordinates whose origin is lo- 
cated at the mass M\ and the angle is measured from 
the separation axis of the binary. The coordinates for 
the source, image and separation vector are denoted by 
(Px,Py) = (p cos ip, psinip), (6 Xl 8 y ) = (r cos (j), r sin 4>) and 
{(-xAy) = (^)O), respectively, where p, r and £ > 0. For 
brevity's sake, cos ip and sin ip are written as C and S, re- 
spectively. The parallel and vertical parts of Eq. (|2|) with 
respect to the separation vector £ are 



r cos 4> ( 1 — 
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r 2 — 2£r cos <fi + £ 2 
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r 2 - 2£r cos <j> + I 2 
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First, let us investigate a case of p = 0, which corre- 
sponds to the source being located behind object 1. Then, 
Eq. (0) becomes 



r cos ( 
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r 2 — 2£r cos <j) + £ 2 



2£r cos cb + f 2 ' 



(8) 



whose right-hand side does not vanish because of v ^ 
for the binary. This leads to sin0 = 0, since the R. H. S. 
of Eq. (0) does vanish. By using the Cartesian coordinates 
(0x,Qy) = (r cos 0,0) and 9 2 —r 2 , Eq. (||) is rewritten as 

el - ml -6 X + £{l -v) = 0. (9) 

All solutions for this equation will be given later by Eqs. 
( |3l| ) and (|3^) as the limit of p — ► 0. In the following, let us 
assume that p ^ 0. As for S, there are two cases, off-axis 
sources (S ^ 0) or sources on the symmetry axis (5 = 0). 

2.1. Off-axis sources 

Here, we consider the case that S does not vanish. In order 
to eliminate the common factor in the L. H. S. of Eqs. (j^j 
and (0), we divide Eq. (§) by Eq. (§), whose R. H. S. does 
not vanish for 5^0. We obtain 

vl 



r 2 -*tr C o S <t> + f = p{c _ Scot(f)y (10) 
Substituting this into r 2 — 2£r cos (f> + £ 2 in Eq. (0) gives 



us 



£ - p(C - S cot <j>) r 2 sm(b - £pSr -£{l-v) sin cb = 0.(11) 

We eliminate r 2 in this equation by using Eq. (|lo|). Hence, 
we obtain 



(12) 
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Ri{4>) = (£ 2 P C-£p 2 C 2 - 
- P S{£ 2 - 2£pC 
-£ P 2 S 2 cos 2 <j>, 
#2(0) = p(C sin (j) — S cos 4>) 

' - pC) sin 2(f) + pS cos 2</> 
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We can show that r cos is a function only of tan < 
Namely, Eq. @ is rewritten as 



r cos ^ 



i?i(tan0) 



i? 2 (tan (j)) 

where we defined 

(£ 2 pC - £p 2 C 2 -vl + P C) tan 2 (j> 
-pS{£ 2 - 2£pC + 1) tan - £p 2 S 2 , 
p(Ctan0- 5) 

pS + 2(£ - pC) tan - pS tan 2 <; 
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Equation ( |15[ ) plays a crucial role in this letter; if we find 
out tan 4>, Eq. ( |l5| ) gives us the value of r cos <f>. The re- 
maining task is deriving an equation for tan <j>. 

Let us substitute Eq. ( p"2| ) into Eq. (pTj|), so that r can 
be eliminated. After lengthy but straightforward compu- 
tations, we obtain an equation for tan</> 



(05 tan 5 



04 tan 
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-I- 03 tan 
+ai tan</> + ao) tan</> = 0, 

where by the frequent use of C 2 + S 2 

a = v£p 3 S 3 , 

01 = /) 2 S 2 + 2VCS 2 -l 2 (2p 2 S 2 -p 4 S 2 ) 

~2£ 3 p 3 CS 2 +£ 4 p 2 S 2 

~v{Up 3 CS 2 - 4^ 2 p 2 S 2 ), 
a 2 = -2p 2 CS-A£p 3 {C 2 S-S 3 ) 

+£ 2 (4p 2 CS - 2 /9 4 CS') + A£ 3 p 3 C 2 S - 2£ i p 2 CS 

+v \£{2pS + 8p 3 C 2 S - 2p 3 S 3 ) 

- m 2 p 2 CS + 2£ 3 pS . 
a 3 = p 2 + £{2p 3 C 3 - \Qp 3 CS 2 ) 
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+£ 2 (-2p 2 C 2 + 2p 2 S 2 + p A ) - 2£ 3 p 3 C 
+v\t(-2pC - 4p 3 C 3 + 6p 3 CS 2 ) 



£ 4 p 2 



+ U 2 p l C z 
+v 2 £ 2 , 

a 4 = -2p 2 CS + 8£p 3 C 2 S - £ 2 {Ap 2 CS + 2p A CS) 
+4£ 3 p 3 C 2 S ~2£ i p 2 CS 

+v \t(2pS - Ap 3 C 2 S + p 3 S 3 ) 

- 2£ 2 p 2 CS + 2£ 3 pS\ , 
a 5 = p 2 C 2 - 2£p 3 C 3 + £ 2 (2p 2 C 2 + p 4 C 2 ) 
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-2£ 3 p 3 C 3 + £ A p 2 C 2 

+v \-£(2pC + p 3 CS 2 ) + 2£ 2 p 2 C 2 - 2£ 3 pC 

We 2 . 



(24) 



It should be noted that all of these coefficients ao, ■ ■ ■ 05 
are not singular, since they are polynomials in £, p, v, C 
and S all of which are finite. In the case of nonvanishing 
p and S, Eq. (0) means that smd> does not vanish and 
consequently neither tan <fi. Hence, Eq. ( |l8| ) is reduced to 
the fifth-order equation for tan 0, 



a,i{ta,YL(j)) % = 0, 



(25) 



As shown by Galois in the 19th century, a fifth-order 
equation cannot be solved in the algebraic manner (e.g. 
van der Waerden 1966). Hence, by solving numerically 
Eq. (|25|), the image position is obtained as (9 X ,6 V ) = 
[r cos cf>, r cos <j> tan (f>) . It is important to consider a re- 
lation of Eq. (||) to the treatment of Witt (1993) in 
which a single complex fifth-order algebraic equation for 
z = 9 X + iOy is obtained: When we use a relation tan <fi = 
—i{z — z)/(z + z), Eq. ( p5| ) can be derived also from the 
complex equation after lengthy manipulations. 

For a source inside the caustic, Eq. (25) has five real so- 
lutions corresponding to five images, while it has three real 
and two imaginary solutions when the source is outside 
the caustic (Witt 1990, Erdl & Schneider 1993, Witt & 
Petters 1993). This criteria can be re-statcd algebraically 
by the use of the discriminant D5 for the fifth-order equa- 
tion (E^), which takes a rather lengthy form in general, 
namely 59 terms (e.g. van der Waerden 1966). It is worth- 
while to mention that all of real solutions for Eq. ( p5| ) must 
exist between —K and K, where K is the larger one be- 
tween 1 and |ao/a 5 | + 1 01/05! + (02/05! + (03/05! + |a 4 /a 5 |. 
(For instance, see the section 66 in van der Waerden 1966). 

In numerical computations, it might be difficult to han- 
dle extremely large r = tan</>. In such a case, we can sep- 
arate |t| € [0,oo) into |t| 6 [0,r c ) and |r| e [t c ,oo). We 
choose t c so that Eq. p^ ) can be easily solved for |r| 6 
[0, t c ). On the other hand, for |r| 6 [r C) 00), instead of Eq. 
(p5|), we can solve an equation for cot <j> 6 (— l/r c , l/r c ) 



i=0 



0. 



(26) 



2.2. Sources on the symmetry axis 

Let us consider the case of vanishing S (C = ±1), namely 
sources on the symmetry axis, for which analytic solutions 
can be obtained: For a binary with two equal masses, ex- 
plicit solutions were found by Schneider and Weifi (1986), 
while no explicit solutions have been given for an arbitrary 
mass ratio until now. Hence, analytic solutions which are 
given below are useful for verification of numerical imple- 
mentations, since numerical solutions for 5* ^ in the 
subsection 2.1 must approach analytic ones as 5* — * 0. 



For S = 0, Eq. (0) implies apparently the following 
three cases 



r = 0, 
sin0 = 0, 

1 - v 
1 w- 



r 2 — 2£r cos (j> + £ 2 



= 0. 



(27) 
(28) 

(29) 



The case of r — should be discarded, since the left-hand 
side of Eq. (0) diverges. 

Next, let us consider the case of sin</> = 0. For this 
purpose, it is convenient to use the Cartesian coordinates 
(e xi 9 y ) = (r cos 0,0) and ((3 x ,f3 y ) = {pC, 0) = (p, 0) for 
p € (— 00,00). By using r 2 = 6 2 , Eq. @ is rewritten as 
the third-order equation 



9 3 x - {£ + p)Bl + (£p - 1)0, + £{l - v) = 0, 

which coincides with Eq. ^ as p — > 0. Equation 
solved explicitly as 




f i cos- 
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where 



P 

q : 



+ p) A + -{£ + p){£p-i) + £{l-v). 



(33) 
(34) 



2 

27 v " ' r/ ' 3' 

Actually, we can show that p < and g 2 + 4p 3 < 0, which 
mean these three solutions exist for any source position. 
They are on-axis solutions, while the remaining two so- 
lutions are off- axis solutions, which are present only for 
sources within the caustics. 

Finally, for the case defined by Eq. ©, Eq. (g) leads 

to 



P 

f 



r 2 - 2£r cos 4> + £ 2 

Replacing r 2 — 2£r cos < 
we obtain 

2 1(1 

which has the positive solution 



in Eq. rt2 



11(1 - v) 
£-p ' 



(35) 

by this equation, 
(36) 

(37) 



if and only if £ > p. Substitution of the solution into Eq. 
I5I) leads to 



(38) 
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A condition | cos <j>\ < 1 is rewritten as 

£ 2 p 4 - 2£(f - 1 + 2v)p 5 + \{1 - I 2 ) 2 + &vf 
-2l{l + l 2 )vp + v 2 l 2 < 0. 



(39) 



Only if this is satisfied for p < £, the two images appear 
at the location given by Eqs. (|37|) and (^3). 



3. Conclusion 

We have carefully reexamined the lens equation for a bi- 
nary system in the polar coordinates. As a consequence, 
we have derived Eq. (|2^) for tan <fi. After solving the equa- 
tion, rcos0 is determined by Eq. (|l5|). Hence, the im- 
age position (8 x .8 y ) — (r cos 0, r cos <fi tan (j>) can be de- 
termined. Our formulation based on the one-dimensional 
equation (^5|) is significantly useful compared with pre- 
vious two-dimensional treatments for which there are no 
well-established numerical methods (Press et al. 1988); the 
new formulation enables us to study the binary lensing 
more precisely with saving time and computer resources. 
For instance, it is effective in rapid and accurate light- 
curve fitting to microlensing events, in particular due to 
star-planet systems. 
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